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Charged particles and neutral kaons in 
photoproduced jets at HERA 
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Charged particles (/i ± ) and K° mesons have been studied in photoproduced events 
containing at least one jet of Et > 8 GeV in a pseudorapidity interval (-0.5, 0.5) in the 
ZEUS laboratory frame. Distributions are presented in terms of transverse momentum, 
pseudorapidity and distance of the particle from the axis of a jet. The properties of 
within the jet are described well using the standard settings of PYTHIA, but the use of 
the multiparton interaction option improves the description outside the jets. A reasonable 
overall description of the K° behaviour is possible with PYTHIA using a reduced value 
of the strangeness suppression parameter. The numbers of and K° within a jet as 
defined above are measured to be 3.25 ±0.02 ±0.28 and 0.431 ±0.013 ±0.088 respectively. 
Fragmentation functions are presented for and K° in photoproduced jets; agreement 
is found with calculations of Binnewies et al. and, at higher momenta, with pp scattering 
and with standard PYTHIA. Fragmentation functions in direct photoproduced events are 
extracted, and at higher momenta give good agreement with data from related processes 
in e + e~ annihilation and deep inelastic ep scattering. 
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1 Introduction 



When a high energy photon interacts with a proton, a QCD-governed hard scattering process 
may take place in which either the photon, or a parton within the photon, interacts with a parton 
in the proton. A signature for such processes is high-.&r jets in the final state. Both HI and 
ZEUS at HERA have reported measurements of various aspects of jet photoproduction [1-5]. In 
general terms, the properties of particles in jets depend on the type of leading parton generated 
by the hard process, and on the process known as fragmentation (or hadronisation) by which 
the colour in the jet is neutralised. The main features of the fragmentation of a given leading 
parton at a given energy are believed to be universal and independent of the type of initiating 
process. It is fragmentation, moreover, rather than leading strange or charm quarks, which 
generates most of the strange particles that occur in jets 0. At the phenomenological level, 
a longitudinal fragmentation function D(z) may be used to describe the observed momentum 
distribution of a chosen type of particle along the jet axis, where z is the fraction of the jet 
momentum taken by the particle. Experimentally, D(z) functions have so far been obtained 
mainly from fits to e + e~ collider data. 

In the present analysis, using data taken in ep collisions with the ZEUS detector in 1994, we 
examine the inclusive properties of charged particles and K° mesons produced in association 
with high-i?T photoproduced jets. This work complements the studies of such particles which 
have been performed in deep inelastic scattering (DIS) by the ZEUS and HI collaborations [7- 
10]. The results are compared with predictions obtained from the PYTHIA Monte Carlo 
generator, which incorporates the Lund string model of particle fragmentation as implemented 
in JETSET [[TTJ. As one possible improvement to the standard description, the suggestion is 



investigated that the particle distributions may be affected by an underlying event structure 
due to multiparton scattering 



Several studies are also made in this analysis to test the universality of fragmentation. In 



deep inelastic scattering j9L hG, 13], it has been found that the behaviour of charged particles 



in the current region of the Breit frame is similar to that found in high energy quark jets in 



e + e~ annihilation; however, there are indications j7|, || 13] that DIS may give jets with a lower 
strangeness content than expected from the standard settings of PYTHIA. DELPHI [15] has 
also reported results with a similar conclusion. We investigate whether this applies also to 
jets produced using incoming photons that are quasi-real, noting that inclusive cross sections 



of K in photoproduction found by HI are broadly consistent with PYTHIA predictions |TB 
Comparison is also made with the fragmentation function calculations of Binnewies et al. [17 



which are based on fits to the particle content of jets at PEP and LEP |Tj|]. To give a more 



immediate test of the universality of fragmentation in different types of process, we select a 
sample of events dominated by the direct photoproduction process (i.e. in which the photon 
interacts in a pointlike manner). These data allow a comparison with e + e~ and DIS results 
without any intermediate model or parameterisation. 

The structure of the paper is as follows. After an account of the apparatus, the event selection 
and the selection of charged tracks, the procedure for reconstructing K° mesons is discussed. 
We describe the Monte Carlo models used in simulating the data, and the more important 
sources of systematic error. Results are then given on the properties of charged particles and 
K° mesons in photoproduced jets. Finally, we present a study of fragmentation functions for 
both charged particles and K° mesons, extracting results for direct photoproduced events in 
order to make a general comparison with other types of experiment. 
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2 Apparatus and running conditions 



The data used in the present analysis were collected by the ZEUS detector at HERA. During 
1994, HERA collided positrons with energy E e = 27.5 GeV with protons of energy E p = 820 
GeV, in 153 circulating bunches. Additional unpaired positron (15) and proton (17) bunches 
enabled monitoring of beam related backgrounds. The data sample used in this analysis cor- 
responds to an integrated luminosity of 2.6 pb _1 . The luminosity was measured by means of 
the positron-proton bremsstrahlung process ep — > cyp, using a lead-scintillator calorimeter at 
Z = —107 mP] which intercepts photons radiated at angles of less than 0.5 mrad with respect 
to the positron beam direction. 

The ZEUS apparatus is described more fully elsewhere |I9| . Of particular importance in the 
present work are the central tracking detector (CTD), the vertex detector (VXD) and the 
uranium-scintillator calorimeter (CAL). 

The CTD [^0[ is a cylindrical drift chamber situated inside a superconducting magnet coil 
which provides a 1.43 T field. It consists of 72 cylindrical layers covering the polar angle region 
15° < 9 < 164° and the radial range 18.2-79.4 cm. The transverse momentum resolution for 
tracks traversing all CTD layers is u(pt)/pt ~ y / (0.005pr) 2 + (0.016) 2 , with px in GeV. The 
VXD [^1| supplied tracking inside the CTD, and consisted of 120 radial cells, each with 12 sense 
wires covering the radial range 10.6-14.3 cm. The vertex position of a typical multiparticle event 
is determined from the tracks to an accuracy of typically ±1 mm in the X, Y plane and ±4 
mm in Z. 

The CAL |22| gives an angular coverage of 99.7% of 47r and is divided into three parts (FCAL, 
BCAL, RCAL), covering the forward (proton direction), central and rear polar angle ranges 
2.6°-36.7°, 36.7°-129.1° and 129.1°-176.2°, respectively. Each part consists of towers which are 
longitudinally subdivided into electromagnetic (EMC) and hadronic (HAC) readout cells. From 
test beam data, energy resolutions of cte/E = 0.18/v^E for electrons and & E /E = 0.35/-\/E 
for hadrons have been obtained (with E in GeV). The calorimeter cells also provide time 
measurements which are used for beam-gas background rejection. 

3 Trigger and event selection 

To identify jets in the event trigger and the subsequent selection of events, a cone algorithm |4[] 
in accordance with the Snowmass Convention p3f was applied to the calorimeter cells. Each 
cell signal was treated as corresponding to a massless particle. A cone radius of 1.0 in R — 
\J (5(f)) 2 + (Sr/) 2 was used, where 8(f), 5r] denote the distances of cells from the centre of the jet 
in azimuth and pseudorapidity. A similar algorithm was used both online and offline. The 
transverse energy of the reconstructed jet is the sum of the measured transverse energies of the 
calorimeter cells included in it, and will be referred to as Ef ec . 

The ZEUS detector uses a three-level trigger system. The first level trigger selected events 
on the basis of a coincidence of a regional or transverse energy sum in the calorimeter and 

^^The ZEUS coordinates form a right-handed system with positive-Z in the proton beam direction and a 
horizontal X-axis. The nominal interaction point is at X = Y = Z = 0. Pseudorapidity r\ is defined as 
— In tan 6/2, where 9 is the polar angle relative to the Z direction. In the present analysis, ij is always defined 
in the laboratory frame, and the actual interaction point of the event is taken into account. 
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a track in the CTD pointing towards the interaction point. At the second level, at least 8 
GeV of transverse energy was demanded, excluding the eight calorimeter towers surrounding 
the forward beam pipe. Beam-gas background is further reduced using the measured times 
of energy deposits and the summed energies in the calorimeter. At the third level, jets were 
identified. The trigger combination used in the present study required at least one jet with 
Ef ec > 6.5 GeV and i] < 2.5, or with Ef ec > 5.5 GeV and r\ < 2.0. Cosmic ray events were 
rejected by means of information from the CTD and the calorimeter. An interaction vertex 
at a position Z > —75 cm, as determined from the CTD tracks, was demanded. A total of 
392k events passed these requirements. The trigger efficiency is close to 100% over the entire 
kinematic range of events used in the present analysis. 

To study the association of K® mesons with jets, it was necessary to obtain a sample of high 
energy jets sufficiently centred in the acceptance of the CTD so as to optimise the acceptance 
for the 7r + and the 7i~ decay products of associated K® mesons. With these considerations in 
mind, events were selected offline with at least one jet having Ef ec > 7 GeV and \r]\ < 0.5. 



Standard ZEUS background rejection criteria were applied [L9| to improve the rejection 
of beam-gas events and cosmic ray events by means of cuts on the primary vertex position, 
the fraction of well-measured tracks, the CAL signal times and the transverse momentum 
imbalance in the event. On the assumption that the outgoing positron is not detected in the 
CAL, Hjb = J2(E ~pz)/2E e was calculated, where the sum is over all calorimeter cells, treating 
each signal as equivalent to a massless particle; i.e. E is the energy deposited in the cell, and pz 
is the value of Ecos8. The quantity yjB is then a measure of y true = E 7 ^ in /2E e , where E 1>in is 
the energy of the incident virtual photon. In the case that a DIS positron is present, a value of 
approximately unity is obtained. Events containing a DIS positron were rejected by requiring 
that (i) no scattered beam positron be identified in the CAL, and (ii) y,js < 0.7. A requirement 
of yjB > 0.15 was also imposed as part of the beam-gas background rejection procedure. A 
total of 34.8k events was obtained at this stage, containing 36.6k jets with Erf ec > 7 GeV and 
\t)\ < 0.5. This selection of events, each containing at least one accepted jet, represents the 
basic event sample for the analysis that follows. 



4 Charged particle selection 

For the study of inclusive charged particle distributions, CTD tracks were accepted if 

(1) they were associated with the primary vertex, 

(2) their pseudorapidity was in the range \r)\ < 1.5, and 

(3) their transverse momentum satisfied pt > 0.5 GeV. 

The pseudorapidity and momentum conditions are chosen to be the same as those to be imposed 
on the mesons. Apart from a small component arising from short-lived particle decays, 
these conditions provide a sample of well-measured tracks which can be identified with charged 
hadrons originating from the main interaction. In the following sections, charged particles will 
be referred to as h^. 
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5 K° reconstruction 



K° mesons decay 50% as K®, which have a 68.6% branching ratio into 7r + n~ ^4|. K® mesons 
were identified by their charged decay mode K® — > 7r + 7r~. This decay mode is easily identifi- 
able given the accurate tracking measurements from the CTD, since the mean decay distance, 
projected on to the r<fi plane, of 2.67pT/mx cm ensures a spatial separation of the decay ver- 
tex from the primary event vertex for a large number of the K® produced in the kinematic 
conditions of the present analysis. 

Kg identification starts by selecting pairs of oppositely charged tracks obtained using the stan- 
dard ZEUS reconstruction algorithms. In order to restrict the track selection to the kinematic 
region where the tracking was best understood, each track was required to satisfy the conditions 
Pt > 150 MeV and pseudorapidity \rj\ < 1.75. More details of the track reconstruction proce- 
dure are given in ref 0. Starting with the track parameters determined at the point where the 
track passes nearest to the beamline, and assuming the tracks in this region to be circular in 
the r0 plane, candidate secondary vertices were then found by first calculating the intersection 
points in the rip plane of all pairs of oppositely charged tracks in an event. Zero or two such 
intersection points are found for a given track pair. The Z-coordinates and the momentum 
components of each track of a pair were then calculated at each of these (r, <p) points. 

Further selections were made on the kinematic quantities listed below. 

(i) Separation in Z position. The separation |AZ| of the two tracks at a candidate secondary 
vertex was required to be less than 3 cm. If a given track pair gave two candidate 
secondary vertices, the one with smaller |AZ| was chosen. 

(ii) Collinearity. The collinearity angle a is defined as the projected angle in the r<p plane 
between the K® momentum and the line joining the primary to the secondary vertex. 
Candidates were accepted if cos a > 0.99. 

(iii) Track impact parameter. The impact parameter e of a track is defined as the distance 
of closest approach between the extrapolated track and the primary event vertex. A 
requirement e > 0.3 cm was imposed on both tracks for each kaon candidate. 

(iv) Photon conversion. The effective mass of the two tracks was evaluated assigning zero 
mass to each track. Track pairs whose effective mass was less than 50 MeV were excluded 
from further consideration. 

(v) A removal. The pre mass hypothesis was applied to each track pair, taking the higher 
momentum particle to be the p (p). As in 0, K® candidates with a pn mass less than 
1.12 GeV were rejected. 

mesons with pt > 0.5 GeV and \rj\ < 1.5 were used in the present analysis. The 7r + 7r~ 
mass distribution for all accepted candidates in this kinematic range is shown in fig. |l|, and 
displays a strong K® signal. After subtracting a background (averaged over windows between 
440-470 and 530-560 MeV) from the signal region of 470-530 MeV, a total of 3154±63 K° was 
evaluated on a background of 388. The mean value of the reconstructed K® mass was 497.0±0.1 



MeV (statistical error) compared with the nominal value of 497.7 MeV [p4|, with a fitted width 
6.3±0.1 MeV. Monte Carlo comparisons show that the decay vertex and momentum of the K® 
are well reconstructed over the full range of momentum and pseudorapidity, without significant 
systematic bias. This background subtraction method was used to obtain the K® signal in 
each bin of all plotted distributions. Further details of the K® reconstruction may be found in 



ref. 251. 
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6 Monte Carlo simulations and data correction 



The experimental data were compared with Monte Carlo generated events obtained using 
PYTHIA 5.7, with JETSET 7.4 for the hadronisation Systematic studies were made 



using HERWIG 5.8 |[26| . PYTHIA is found to give a reasonable description of jet profiles in 
hard photoproduction in the kinematic region of the present study p7| |28fl . PYTHIA events 
were generated using a pr mm value of 2.5 GeV as the cut-off for the hard subprocess; the results 
were insensitive to this choice. The event generation was followed by a full simulation of the 
detector and trigger response in ZEUS by means of GEANT 3.13 |2l| . The direct and resolved 
contributions are combined in proportion to their generated cross sections. The GRV-LHO 
parton densities for the photon |3D| and MRSA for the proton |JT]] were used in running the 



standard version of PYTHIA. Variations on the standard calculation were made by reweighting 
the results to use different photon parton densities, and also by using an option which uses 
multiparton interactions (MI) as a model for an underlying event accompanying the main hard 
QCD subprocess (with a Pt mm value of 1.4 GeV for the secondary interactions). In the LO 
QCD model employed here, multiparton scattering may accompany the resolved processes, in 
which the photon acts as a source of partons, but cannot accompany the direct process, in 
which the photon interacts as a pointlike entity with no partonic substructure. 

Within string fragmentation models such as PYTHIA, the most important quantity governing 
the number of strange particles that appear in the fragmentation is the so-called "strange- 
ness suppression parameter" P s /P u . This is assigned a default value of 0.3 on the basis of 



measurements at PETRA and PEP |32|], which measured a variety of strange/non-strange 
particle ratios (e.g. K° to mesons) in the products of e + e~ annihilation. In previous studies 
of neutral kaons in DIS at HERA |7], [J it was found that decreasing P s /P u from its standard 
value improved the agreement between the Monte Carlo and the data. We have therefore 
generated Monte Carlo event samples also with the value of P s /P u decreased from 0.3 to 0.2. 

The Monte Carlo events were used to determine correction factors for the data. These were 
calculated separately for each bin of any given plot, so that for each plotted quantity the 
corrections are suitably averaged over the other physical quantities. For K° calculations, the 
following factors were taken into account: 

(i) K® reconstruction efficiency e^o ■ Given a set of reconstructed and accepted Monte Carlo 
events, this was the ratio of (a) the number of reconstructed — > 7t + tt~ entries in a given 
bin, background subtracted, to (b) the corresponding number of generated — > 7r + 7r~ 
occurring in this bin. All the other event selection criteria are applied in the normal way. 

(ii) Reconstruction level/hadron level correction factor C . This factor makes use of jets re- 
constructed from the four-vectors of the primary final state particles (charged or un- 
charged) in the generated events, referred to as "hadron jets" . For events with a generated 

— > 7t + tt~ in a given bin, C is defined as the ratio of (a) the number of events which 
satisfy the trigger and acceptance conditions and have a reconstructed calorimeter jet 
with E^ ec > 7 GeV and \rj\ < 0.5, to (b) the number of events having 0.2 < y true < 0.85 
and a hadron jet with cone radius 1.0, Ej- > 8 GeV and \r]\ < 0.5. In this way we correct 
to a defined set of kinematic conditions at the final state level. The parameters in (b) are 
chosen to provide a good correspondence with the cuts at the detector level, and thereby 
minimise the corrections. The y true range corresponds to a ^p centre of mass energy range 
of 134-277 GeV. 
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Since the shapes of the distributions are similar in the data and Monte Carlo (after recon- 
struction), the calculated correction factors will suitably take into account all effects due to 
reconstruction efficiency, event selection efficiency and bin-to-bin migration. These aspects 
have been checked in our previous studies of inclusive jet production Q. The K® reconstruc- 
tion efficiency is dominated by the geometric effects of the track cuts and the CTD tracking 
efficiency, with little sensitivity to the properties of the jet. It is similar for direct and resolved 
events. In the region of px and r\ used here, e K o has a plateau of approximately 0.35 in the 
middle of the accepted range of either quantity, falling to 0.25 at the ends. C takes values 
typically between 0.6 and 0.8. 

Corrected numbers of K mesons in bins Av of any given variable v are evaluated according to 
the formula 

dN (corrected) N(detected) 

lb = Ave tot B ' ( } 

where e to t = £r° C, and B is the total branching ratio K° — > K® — > 7r + 7r~. 

For charged particles, a similar relationship e toi = e^± C is used, where the charged particle 
reconstruction efficiency e%± is calculated as the ratio of reconstructed and accepted tracks to 
the number of generated charged particles in the same bin. C is then calculated analogously 
to the K° case, and B is set to unity. Here and in the calculation of C, generated charged 
particles are taken to be those in the final state with lifetimes greater than 10 -8 s, together with 
charged decay products of primary hadrons with a shorter lifetime than this, but excluding the 
products of K® and A, A decays. 

All the particle distributions are normalised to the number of jets. Keeping in mind that each 
event in the basic event sample has at least one accepted jet, we evaluate (a) the total number 
of charged particles or kaons in each bin of a given distribution, and (b) the total number of 
accepted jets in the basic event sample. The ratio of (a) to (b) is then plotted to give numbers 
of particles per jet. To obtain corrected distributions, it is thus necessary to correct the total 
number of measured jets Nj ets , in the defined kinematic range, to the total number of above- 
defined hadron jets. To achieve this, the correction formula (1) is again employed, using an 
efficiency factor etot — C together with B — 1, and removing from the definition of C the 
conditions on kaons or charged particles. 



7 Systematic uncertainties 

The stability of the results against reasonable variations of the selection criteria was studied in 
order to estimate the systematic uncertainties. The most important effects were estimated as 
follows: 

(1) To study the sensitivity to the fragmentation scheme in the Monte Carlo, HERWIG was 
used instead of PYTHIA to evaluate the correction factors. 

(2) An uncertainty exists in the matching of the calorimeter energy scale in the Monte Carlo 
events to that of the data. To estimate this, the Monte Carlo energy scale was varied by 
±5%. 

(3) To estimate tracking uncertainties, the effect of including or not including the VXD in 
the track reconstruction was investigated. 
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(4) The effect of varying the accepted i/jb range to 0.2 < jjjb < 0.8 was evaluated. 

(5) K° definition: 

(i) the cut on \AZ\ was varied between 2 and 4 cm; 

(ii) the cut on cos a was varied by ±0.005; 

(iii) the cut on the track impact parameter |e| was varied in the range 0.27-0.33 cm; 

(iv) the upper mass value for A rejection was varied between 1.117 and 1.123 GeV. 

The effects on the results of most of the parameter variations listed were normally found to 
be small, at the level of up to a few percent. The individual contributions were combined 
in quadrature to give the total systematic error on a given point. The systematic errors on 
the correction factors are dominated by the effects of (1), (2) and (3). The comparison of 
uncorrected data with PYTHIA predictions is affected by (2) and (3). 

8 Results 

8.1 Charged particles 

In fig. ^ we present efficiency and acceptance corrected distributions of dN^h^/dp? and 
dN^h^/drj in prih^) and r)(h ) respectively. As defined in section 6, the corrected distri- 
butions are given for particles in events containing a final state hadron jet of _&r(jet) > 8 
GeV with |^(jet)| < 0.5. Here and throughout, all results are normalised as numbers of par- 
ticles per jet satisfying the stated definitions. The prih^) distributions are averaged over 
— 1.5 < rjiji^) < 1.5; the ^(h^ distributions are averaged over Pt(^ ± ) > 0.5 GeV. At this stage 
no explicit association of the h with the jets is made. The statistical errors are small. 

Good agreement is seen between the data and the standard PYTHIA predictions in the shape 
and magnitude of the px distribution. In the t] distribution, the overall agreement is still 
good, but it is apparent that standard PYTHIA tends to undershoot the data at higher rj 
values. A similar trend has been observed in other photoproduction studies at HERA 0, f|, |2"7|] , 
leading to suggestions that the discrepancy can be attributed at least in part to multiparton 
interactions ||12|| . The MI option is seen to give an improved description of the data at high 
rj; however, it overestimates the exponential slope of the pt distribution. A reduction in the 
Ps/Pu value has a negligible effect on the present distributions (not shown). Little difference is 
seen if different photon parton densities are used. 

We now examine in more detail the association of charged particles with jets, given the presence 
of at least one jet in each event in the selected data sample. Fig. |3|(a) shows the uncorrected 
distribution in the difference A0 in azimuth between an h and the axis of any jet in the event 
satisfying Erf ec > 7 GeV and |?7(jet)| < 0.5. The peak around |A0| = 7r indicates the presence 
of jets opposite in azimuth to the observed particle. For h -jet pairs with |A0| < n/2, the 
distance Ar] in pseudorapidity between the and the jet is plotted in fig. ||(b). In comparing 
the data to the Monte Carlo predictions, an overall systematic uncertainty of ±5% should be 
allowed. 

In discriminating between the various models, a particularly useful variable was found to be 
R = J (A</>) 2 + (Ar/) 2 , i.e. the distance in (rj, <p) between a given particle and the jet axis. The 
uncorrected R distribution is presented in fig. [|(c). A prominent peak at small R demonstrates 
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that the are being produced in association with jets; there is also a peak at R ~ 3 due 
to the likely existence of another jet opposite in to the first. In the following discussion 
we concentrate mainly on the region R < 2.5. Over the whole range, the main contributions 
come from resolved processes. The predictions of standard PYTHIA lie slightly above the data 
at low R, but are low in the inter-jet region 1 < R < 2.5. To establish whether the latter 
discrepancy is associated with the excess of h over expectations seen at high ?7(/i ± ) a similar 
plot was made with the h acceptance changed to —1.5 < T)(h ) < 0.5 (not shown). This 
reduced the discrepancy for 1 < R < 2.5 slightly but did not eliminate it. Variations of the 
P s /P u parameter and the photon parton density again have little effect on the distributions. 
The effect of using the MI option of PYTHIA is more pronounced; in the inter-jet region, the 
MI option overcompensates for the previous shortfall of standard PYTHIA events, while at 
small R the jet appears broadened. Overall, the MI option provides an improved description 
of the data. 

The mean corrected multiplicity of charged particles with p^(/i ± ) > 0.5 GeV per jet, integrated 
over R < 1, is shown in Table [l] compared with different results from PYTHIA. All the models 
are consistent with the experimental result, within errors. 

8.2 K° mesons 

In the previous section it was seen that the features of charged particles in photoproduced 
jets are generally well described using standard PYTHIA, with some further improvement 
obtainable using the MI option. Turning now to K° mesons, we consider the same distributions 
as for the h . As before, at least one reconstructed calorimeter jet with Ef ec > 7 GeV 
and |^(jet)| < 0.5 was demanded in an event, and the corrected distributions are for events 
containing a final-state hadron jet with .Bp (jet) > 8 GeV and |?y(jet)| < 0.5. The distributions 
are normalised as numbers of K® or K° per jet for uncorrected or corrected data respectively. 
Corrected dN(K°)/dp£ and dN(K°)/dr} distribut ions are shown in fig. [|as functions of px(K°) 
and i](K°). The pt{K°) distributions are averaged over —1.5 < rj(K°) < 1.5; the 7)(K°) 
distributions are averaged over pt(K°) > 0.5 GeV. No attempt is yet made to associate the 
Kg with jets. In the larger error bars the systematic and statistical errors are summed in 
quadrature. Following the studies mentioned |7|, [J, we also show results with the value of 
P s /P u decreased from 0.3 to 0.2. 

A reasonable overall agreement is seen between the data and the various PYTHIA plots in the 
shape and magnitude of the pt spectrum, although the MI option tends to overestimate the 
slope of the plot. For rj(K ) < 0.5, the data of fig. |](b) show an approximate agreement with 
the standard PYTHIA which is improved by using a reduced value of P s /P u . At higher rj(K°), 
however, a tendency for the PYTHIA values to be low compared with the data is worsened 
when P s /P u is reduced. The agreement in this region is improved with the use of the MI option. 

Fig. |5](a) shows the uncorrected distribution in the difference A0 in azimuth between a K® 
and any jet in the event satisfying E^ ec > 7 GeV and |r/(jet)| < 0.5. For K^-jet pairs with 
\A<fi\ < 7r/2, the distance A77 in pseudorapidity between the and the jet is plotted in fig. 
0(b). A common overall systematic uncertainty of ±7% should be allowed. The data are 
reasonably well fitted by the Monte Carlo outside the peak at A<ft « and in the wings of 
the T)(K°) distribution. In both distributions, standard PYTHIA overestimates the numbers of 
kaons near the axis of a jet. 
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The uncorrected distribution in R is plotted in fig. |5|(c) for all Kg. Here it is clear that standard 
PYTHIA overestimates the data in the jet core, i.e. for R < 0.4, while underestimating it in 
the interjet region. To gain understanding of the possible reasons for the discrepancy in the jet 
core, a variety of investigations were made. It was found that the discrepancy persists when 
the number of charged tracks in the region of the jet (i.e. with R < 1) is selected to be small: 
even with just one or two charged particles present in addition to the K®, a similar effect is 
seen. In the events used in this analysis, the mean number of additional charged particles in 
the jet is approximately 3. It is therefore difficult to attribute the effect to a poorly understood 
problem with the charged particle tracking. 

The effects of varying the photon parton densities were investigated and found to be small: 



the use of the GRV-LO ||, ACFGP @, LAC1 §§ or GS-HO @ parton sets altered the 
predictions by less than 5% for R < 1. No attempt was made to vary the parton densities in the 
proton, since these are well defined in the present kinematic range from other measurements. 
Removing the small contribution from charm-containing jets in the generated events made 
only a slight difference, giving little scope for remedying the problem by remodelling the charm 
simulation. If the normalisation is performed to luminosity rather than to numbers of jets, 
the conclusions are likewise unchanged. A similar although less marked discrepancy at low 
R was found using predictions from the HERWIG Monte Carlo. The discrepancy outside the 
jet is again found to be reduced slightly but not eliminated if the acceptance is reduced to 
— 1.5 < r](Kg) < 0.5. The MI option gives an improved description in the inter-jet region, but 
is worse elsewhere. 

In about 7% of the kaonic events, two were found. The distance R in (jj, </>) between them 
was plotted, giving a distribution which was fairly flat up to R ~ 3.5 with a small enhancement 
at low R. Standard PYTHIA tends to overestimate the data by approximately the same amount 
as at R < 0.4 in fig. ||(c), but over the entire R range. The conclusion is that many of the K° 
pairs are not strongly correlated. 

The characteristics of the K® in jets can be further investigated by distinguishing between direct 
and resolved photoproduction processes, as defined at leading order in QCD. A subsample of 
events was chosen in which at least two jets were found, one satisfying the above experimental 
jet definition and a second having Ef ec > 7 GeV as before, and //(jet) < 2.5. The two highest 
Et jets satisfying these conditions were used to estimate the fraction of the photon energy which 
takes part in the hard subprocess. This estimate is given in terms of the measured parameters 
of the two jets as i° BS = (E^\ c e~ Vl + E^ e 2 c e~ m ) /2E e y JB . For x° BS < 0.75 the event samples are 
dominated by the resolved process, above this value the direct process is more important 



In figs. H(a), 0(a), R distributions for the resolved-enhanced and direct-enhanced data samples 
(x° BS < 0.75, x° BS > 0.75) are shown, and compared with PYTHIA predictions. A small 
admixture of generated direct events is found in the PYTHIA samples with x° BS < 0.75, and a 
relatively larger admixture of generated resolved events in the samples with x° BS > 0.75. In the 
inter-jet region, a deficit in the standard PYTHIA predictions is seen in the resolved-enhanced 
sample but not in the direct-enhanced sample, suggesting the possibility of an association with 
multiparton interactions. In the jet core, the resolved-enhanced sample shows a much larger 
discrepancy between the data and standard PYTHIA than does the direct-enhanced sample. 
The statistical errors on the PYTHIA histograms are similar to those on the data. 

In fig. |6](b), we consider the effects of the reduced P s / P u value and the MI option in PYTHIA on 
the resolved-enhanced event sample. The use of the MI option is manifestly advantageous, and is 
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in fact essential in order to obtain reasonable agreement with the data when using P s /P u = 0.2. 
A P s /P u value between 0.2 and 0.3 appears preferred. While not perfect, the agreement using 
P s /P u = 0.2 together with the MI option represents a significant improvement on the standard 
version of PYTHIA. Similar conclusions emerge from the direct-enhanced sample (fig. 0) in 
which, overall, the reduced P s /P u value significantly improves the match with the data in the 
R distribution. The MI option, affecting only the resolved contribution to the histograms, has 
a positive effect. A P s /P u value intermediate between 0.2 and 0.3 might represent a further 
improvement. 

The mean corrected multiplicity of K° with pt{K°) > 0.5 GeV per jet, integrated over R < 1, 
is shown in Table [I]. The full event sample containing one or more jets per event is used. 
A comparison with the different results from PYTHIA confirms the preference for a lower 
strangeness suppression parameter in describing the overall numbers of K° within the jets. 

8.3 Fragmentation functions 

In the previous sections, angular distributions of charged particles and K® relative to jets were 
studied, averaged over the particle momenta. The momentum distributions of particles within 
a jet may be described by fragmentation functions D(z), where z is a measure of the fraction 
of the jet momentum taken by the particle. D(z) is defined as (1/Nj ets ) dN(X)/dz, where X 
denotes a charged particle or K°. A particle is defined here as being associated with a given 
jet if it satisfies the criterion R < 1. 

The variable z may be defined in several ways. In the usage of CDF |36|| , the longitudinal 
component of the particle momentum along the jet axis is scaled by the jet energy according 
to a formula which may be written as z^ = p(X) ■ n(jet)/£'(jet). Here p denotes the mo- 
mentum 3-vector of a particle, and n(jet) is a unit vector along the jet axis. An alternative 
approach is to take the energy ratio ze = E(X)/E(jet). This is analogous to the definition 
z = E(X) / E(beam) normally used in e + e~ collider experiments, where E(beam) corresponds 
to the maximum possible energy E(max) of a leading parton. One notes that in the latter 
measurements, since they are carried out in the e + e~ centre of mass frame, the calculated 
value of z does not depend on the actual production angle of the final state particle, and so 
is independent of many details of the hadronisation mechanism and higher order effects. No 
explicit identification of jets is required, or association of particles with jets. 

Here we present our results using both definitions of z. The numbers of h , K° and jets were 
corrected as described in section 6. A small systematic overestimate of the z scale at low values 
exists due to losses of low momentum particles from the reconstructed jet. The loss of energy 
in dead material in front of the CAL also has the effect of increasing the measured z value. A 
correction is applied to remove these effects. Results are plotted at the mean z value of the 
events in a series of chosen z intervals. Values of ze are typically 10-20% larger on average 
than those of zl. 

Figs. |8|(a) and [](a) show the resulting D(zl) distributions for and K° respectively, with 
pion masses assigned to the h ± . Also plotted are values of the corresponding distribution 
obtained using PYTHIA. The distributions are averaged over f/(jet) and £r(jet) within the 
same selected ranges as above. For zl < 0.05 the effects of the 0.5 GeV cut-off on pr in the ^ 
and K° data become important. Good agreement is seen between the data and the standard 
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PYTHIA distributions, in the K° case for zl values above 0.15; in this range it is not possible 
to discriminate between P s /P u values of 0.3 and 0.2. The discrepancy noted in the K® R plots 
appears to be concentrated at zl < 0.15, where the bulk of the statistics lie, apart from the two 
highest Zl points with large errors. The MI is in agreement with the h data at low momenta 
now that the association of particles with jets has been made. Also shown in fig. |8](a) are 
D(zl) values from CDF p6| , for in jets of energy ^40-100 GeV produced in pp scattering. 
Although the CDF energies are considerably higher than those of the present measurements, 
and the mixture of leading final state partons is also likely to be different, there is a clear 
similarity of the fragmentation functions in the range 0.15 < zl < 0.7. 

For a comparison with e + e~ data, results from the present measurements are first compared 
with calculations from the next-to-leading order phenomenological fits of Binnewies et al [|I~7f 



(figs. ||(b), |9|(b)). These authors have made use of a variety of e + e~ data sets at different 
energies to extract separate fragmentation functions for the different types of primary parton 
in the range 0.1 < z < 0.8. The approach is based on the different quark couplings to photons 
and to Z bosons, and on topological properties of events with hard gluon emission. A definition 
z = E(X) / E(bea.m) is used; we therefore adopt here the corresponding definition z = ze = 
E(X)/E(jet). Since the primary partons in the photoproduced jets are unidentified, and we 
currently lack a full NLO Monte Carlo simulation of the 'jp process, a shaded region is drawn 
to cover the spread of the calculated D(z) values for the fragmentation of gluons and different 
types of quark (i.e. u, d, s, c). The calculation has been performed using a QCD scale of 
8 GeV, corresponding approximately to the present jet E? values. Details of which parton 
fragmentation functions contribute to the upper and lower bounds of the shaded regions are 
given in the figure captions. 

Good agreement with the calculations is found, to within their uncertainty as applied to the 
present data. A few remarks need to be held in mind: 

(i) Ref. [O includes charged decay products of long-lived primaries such as K®, A in the 



definition of a charged particle. This gives ~ 10% more charged hadrons compared with 
the present method. 

(ii) As stated above, a jet definition, and the requirement that the given particle be associated 
with a jet, are required in our present approach, but not in the e + e~ based analyses. This 
tends to lower the distributions in particular at lower ze values. 

(iii) The experimental ze (or z£) value is corrected relative to the energy of hadron jets. A 
correction to the energy of the the final-state leading parton in the hard process has not 
been attempted here. 

More precise comparisons are available using the direct photoproduction process, since a close 
similarity is expected in general between event properties in direct hard photoproduction, e + e~ 
annihilation, and DIS at low x, where x is the Bjorken variable. At energies comparable to those 
of the present data, e + e~ annihilation is governed by a single electromagnetic vertex, summed 
over the different types of quark, while the LO direct photoproduction process is dominated by 
photon gluon fusion, whose diagrams feature a similar electromagnetic vertex accompanied by 
a quark-gluon vertex which does not depend on the quark type. In a similar way, in DIS events 
at low x, the virtual photon couples to quarks that come predominantly from the sea, having 
evolved through gqq vertices which are likewise flavour independent. 

We therefore make a comparison between fragmentation functions obtained from the direct 
photoproduction process in ZEUS, and values obtained from e + e~ and DIS measurements. 
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Figs. [10], [TT| show results obtained from the present direct-enhanced event samples, corrected 
by means of PYTHIA on a bin-by-bin basis to remove the resolved component so as to be 
equivalent to results from "pure direct" event samples. These are compared with predictions 
from PYTHIA, and with e + e~ measurements at similar centre of mass energies to those of the 
present hard subprocess, i.e. twice the transverse energy of the selected jets. A factor of 0.5 is 
applied to the published e + e~ data to allow for the dominant qq final state. 

For the h data, a further comparison is made to fragmentation functions from ZEUS calculated 



in the current region of the Breit frame in deep inelastic scattering [fH] |. Here, data points are 
used covering a photon virtuality range of 160 < Q 2 < 320 GeV 2 with 0.0024 < x < 0.01. 
The corresponding jet energy is Q/2, i.e. 6.3-8.9 GeV, roughly equivalent to the jets of the 
present events. It should be noted that the TASSO data |37j include all charged products from 
K® and A decays, while the ZEUS DIS and ZEUS 1994 photoproduction data exclude all such 
products; these effects contribute at the level of ±10%. The h ± data compare well also with 
pp data from the ISR (not shown) although a different parton mixture is expected in the 
final state here. For the inclusive K° data of p9| , |40"f at different fixed centre of mass energies, 



the authors quote scaling cross sections which have been converted to fragmentation functions 
to compare with the present data. 

At low ze, the distributions may be affected by differences between the colour flows in the 
different types of event, as well as by the cut-off imposed at 0.5 GeV in p T in the present 
h and K° data, making comparisons between the different data sets difficult. However for 
ze > 0.1 for h , and ze > 0.15 for K°, the present results are in good agreement with the 
standard PYTHIA predictions, and the D(ze) values from the different measurements are also 
in good agreement with each other. This well illustrates the universality that is believed to be 
a property of the quark fragmentation process. 



9 Summary and conclusions 

We have studied the properties of charged particles (h^) and K° mesons in photoproduced 
events in the ZEUS detector at HERA. The K° mesons were studied in the K® — > 7r + 7r~ 
decay mode. In each event at least one reconstructed jet was required in the calorimeter with 
measured Ef ec > 7 GeV, centrally produced in the laboratory frame. The distributions of the 
h^ 1 and K° in these events were studied as a function of a number of kinematic variables; the 
distance R in (rj, 0) of the particle from the axis of a jet displayed information which was not 
so clearly evident from the other distributions. 

Correction factors were applied to evaluate the numbers of and K° per jet at the final state 
hadron level, with E'j'(jet) > 8 GeV and | ?7 (jet ) | < 0.5, for events in the 7p centre of mass 
energy range 134 < W < 277 GeV. Corrected distributions in transverse momentum p? and 
pseudorapidity 77 are given. The corrected numbers of h and K° within a jet are evaluated, 
with a particle defined as being inside a jet if it is within unit radius of the jet axis in (rj, <ft). 
Fragmentation functions D(z) for and K° in photoproduced jets have also been determined. 
In comparing the present results with those from theory and from other experiments, two 
definitions of z are employed, in which the longitudinal momentum component of the particle 
along the jet axis is scaled to the jet energy (z L ), or the particle energy is similarly scaled (ze)- 
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The distribution of h within photoproduced jets is found to be fairly well described by the 
standard version of PYTHIA, whereas that of K° mesons is not. Outside the jets, more 
and K° are found than predicted; the latter situation can be to some extent remedied 
by using a version of PYTHIA which includes a simulation of multiparton interactions in 
resolved events. Taken overall, the numbers of K° within the jets correspond to a reduced 
value of the strangeness suppression parameter P s /P u in PYTHIA. When the data are divided 
into samples enriched in the resolved and direct photoproduction processes respectively, this 
statement remains true in either case. However the effect is concentrated at z values below 
0.15; at higher values the default parameters of PYTHIA give a satisfactory description of 
the data. This suggests a need for further study of fragmentation at low z, where large numbers 
of particles are found but where their association with jets may be less well defined than at 
high z. 

The fragmentation functions determined using inclusive photoproduced jets are found to be in 
agreement with calculations from Binnewies et al. to within the uncertainty of the calculation 
as applied to the present data. A close similarity is seen between the present h^ 1 fragmentation 
functions in the range 0.15 < zl < 0.7 and those observed in pp scattering. Fragmentation func- 
tions have also been extracted for and K° in direct photoproduced jets, and are compared 
with corresponding data from e + e~ annihilation and from deep inelastic scattering. Agreement 
is good for z E > 0.1 and ze > 0.15 for Yi^ and K° respectively. This, together with the agree- 
ment found in this region with PYTHIA, represents a confirmation of the idea of a universally 
valid description of parton fragmentation. 
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DATA 




PYTHIA Monte Carlo 






(ZEUS 1994) 


Standard 


MI 


Standard 


MI 


Ps/Pu 




0.3 


0.3 


0.2 


0.2 




3.25±0.02±0.28 


3.31±0.03 


3.44±0.06 


3.37±0.04 


3.46±0.04 


K° 


0.431±0.013±0.088 


0.513±0.013 


0.560±0.025 


0.416±0.014 


0.426±0.013 



Table 1: Corrected multiplicities of charged particles (h ± ) and neutral K° mesons (pr > 0.5 GeV) per 
photoproduced jet (£r(jet) > 8 GeV, |7y(jet)| < 0.5), cone radius = 1, for 134 < W < 277 GeV. Comparison 
is made between ZEUS data and PYTHIA Monte Carlo predictions with strangeness suppression parameter 
(P s /P u ) = 0.3 (standard value) and 0.2, without and with the multiparton interactions (MI) option. The first 
error quoted is in each case statistical, the second error on the data is systematic. 
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Figure 1: Reconstructed ir + n mass distribution at selected secondary vertices, compared with Monte Carlo 
prediction with equal normalisation. 
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Figure 2: Corrected distributions of charged particles per jet in events containing a hadron jet with i^x(jet) > 8 
GeV and 1 77 (j et ) | < 0.5. (a) (l/N ]ets )dN(h ± )/dp^ as a function of p T (h ± ) for |j?(ft ± )| < 1.5; (b) as a function 
of r](h ) for pr(h ) > 0.5 GeV. Error bars are statistical and systematic combined in quadrature, and are 
dominated by the systematic. The PYTHIA MI Monte Carlo uses a multiparton interaction model. 
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Figure 3: Distance of charged particles (h 1 * 1 ) from jet (a) in </>, (b) in 77, for \A(j>\ < n/2, (c) in {r},4>), with 
comparison to various Monte Carlo predictions as indicated. The entries are uncorrected numbers of charged 
particles (p T > 0.5 GeV, \r]\ < 1.5) per measured jet (E^ ec > 7 GcV, < 0.5) per unit interval of the plotted 
quantity. The histograms correspond to PYTHIA predictions using the GRV-LHO photon structure and default 
strangeness suppression. The direct contribution to the standard PYTHIA calculation is also shown separately. 
A common systematic error of 5% should be allowed on the data points in comparing with the Monte Carlo. 
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Figure 4: Corrected distributions of K° mesons per jet in events containing a hadron jet with ^(jet) > 8 
GeV and |r/(jet)| < 0.5. (a) (l/N jets ) dN(K a )/dp* as a function of p T (K°) for \r){K°)\ < 1.5; (b) 
(1/Nj e ts) dN(K°)/dri as a function of f](K°) for pt(K°) > 0.5 GeV. Inner and outer error bars are statis- 
tical errors and statistical combined with systematic in quadrature. 
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Figure 5: Distance of from jet (a) in <fi, (b) in 77, for \A(j)\ < n/2, (c) in (77, cf)), with comparison to various 
Monte Carlo predictions as indicated. The entries are uncorrected numbers of K® (px > 0.5 GeV, \r/\ < 1.5.) per 
measured jet (Erf ec > 7 GeV, < 0.5) per unit interval of the plotted quantity. The histograms correspond 
to PYTHIA predictions using the GRV-LHO photon structure and default strangeness suppression. The direct 
contribution to the standard PYTHIA calculation is also shown separately. A common systematic error of 7% 
should be allowed on the data points in comparing with the Monte Carlo. 
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Figure 6: Uncorrected R distributions for resolved-enhanced events (:r° BS < 0.75) compared (a) with standard 
version of PYTHIA (P s /P u = 0.3), showing also the direct and resolved contributions separately, (b) with 
PYTHIA using P s /P u — 0.2, MI option, and both variants together. Kinematic definitions are as in fig. [| 
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Figure 7: Uncorrected R distributions for direct-enhanced events (x° BS > 0.75) compared (a) with standard 
version of PYTHIA (P s /P u = 0.3), showing also the direct and resolved contributions separately, (b) with 
PYTHIA using P s /P u = 0.2, MI option, and both variants together. Kinematic definitions are as in fig. [| 
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Figure 8: Corrected fragmentation functions D(z) for charged particles [h ) with pt(/i ± ) > 0.5 GeV, in 
hadron jets with St (jet) > 8 GeV and 1 77 ( j et ) | < 0.5; systematic errors dominate. In (a) comparisons are made 
with variations on the standard version of PYTHIA, and with CDF results in the approximate jet energy range 
40-100 GeV; the dominant systematic errors are shown. In (b) the data are compared with a shaded region 
denoting upper and lower bounds from the phcnomenological fit of IlTj] . The upper bound corresponds to the 
fitted D(z) values for s, d, u quarks, respectively, in the intervals z = 0.1-0.25, 0.25-0.45, 0.45-0.8. The lower 
bound corresponds similarly to u, c quarks in the intervals z — 0.1-0.35, 0.35-0.8, respectively. 
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Figure 9: Corrected fragmentation functions D(z) for K°. Details as for previous figure. The upper bound of 
the shaded region in (b) corresponds to fitted D(z) values from |1?J for d(= s), c, u quarks in the intervals z = 
0.1-0.2, 0.2-0.45, 0.45-0.8, respectively. The lower bound corresponds to <?, d(— s) in the intervals z — 0.1-0.45, 
0.45-0.8, respectively. Thick error bars denote dominant statistical errors, thin denote dominant systematic 
errors. 
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Figure 10: Fragmentation functions D(ze) for h^, obtained from direct-enhanced events and corrected to 
"pure direct" values. Also plotted are (i) standard PYTHIA predictions for the LO direct process, (ii) data 
from e + e~ annihilation at similar jet energies with centre of mass energies as indicated |37j, halved to take 
account of the dominant qq final state, and (iii) data from ZEUS measured in deep inelastic scattering [ [To| with 
Breit frame jet energies in the range 6.3-8.9 GeV. The variable ze is defined as i?(ft. ± )/i?(jet) for ZEUS 1994 
photoproduction data and Ze — E(h ± )/E(m&x) for the other experimental points. Statistical and total errors 
on the ZEUS 1994 data are shown as for previous figures; errors are total errors for the TASSO and ZEUS DIS 
points. 
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Figure 11: Fragmentation functions D(ze) for K°, obtained from direct-enhanced events and corrected to 
"pure direct" values; ze definitions correspond to those for h . Also plotted are standard PYTHIA predictions 
for the LO direct process and data from previous experiments at similar jet energies, with centre of mass energies 
as indicated ]39|, p0| . Errors on the TASSO and HRS points are total errors. 
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